Supersymmetric superheavy dark matter 
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We propose the lightest supersymmetric particle (LSP) as a well-suited candidate for superheavy 
dark matter (SHDM). Various production mechanisms at the end of inflation can produce SHDM 
with the correct abundance, Hlsp'i^ ~ 0.1, if its mass is sufficiently high. In particular, gravita- 
tional production requires that the mass rriLSP of the LSP is above 3 x 10^^ GeV. Weak interactions 
remain perturbative despite the large mass hierarchy, mLSP 3> mz, because of the special decou- 
pling properties of supersymmetry. As a result the model is predictive and we discuss the relevant 
cosmological processes for the case of a superheavy neutralino within this scheme. 

PACS numbers: 95.35.-|-d, 12.60.Jv, 14.80.Ly 



I. INTRODUCTION 

Uncovering the nature of dark matter (DM) is one of 
the most pressing problems of current research in par- 
ticle physics and cosmology. A wealth of observational 
data suggests that a viable DM candidate has to be non- 
baryonic and should be non-relativistic, at least from the 
time of matter-radiation equilibrium on [1] . The various 
particles proposed as cold dark matter (CDM) candidates 
X can be divided in two main sub-categories: Thermal 
relics were at least once during the history of the Universe 
in chemical equilibrium with the thermal plasma, while 
non-thermal relics have either sufficiently small interac- 
tions or a high enough mass mx to be never produced 
efficiently by processes like, e.g., e^e+ XX. 

The present relic abundance fix of a thermal relic 
scales approximately as fix oc l/uann with its annihi- 
lation cross section (lann- Moreover, unitarity of the 
S matrix restricts annihilations into the l.th partial- 
wave of particles with relative velocity v^d as crinn < 
{21 + l)4:TT/{vreiM'^) 0]. Since for non-relativistic point- 
particles higher partial-waves are suppressed, the ob- 
served value [3j] flcDuh^ = 0.11 of the DM abundance 
constrains the mass of any thermal relic as m < 100 TeV. 
On the other hand, the requirement that the DM is cold 
translates for thermal relics into a lower mass limit of the 
order lOkeV. Thus the mass of thermal relics should lie 
in the 10 keV - 100 TeV range. 

Supersymmetry (SUSY) provides with the lightest 
neutralino one of the most attractive candidates for ther- 
mal DM, see Ref. [Jj for a review. Low-energy SUSY 
models are a natural extension of the standard model 
(SM), offering a solution to the hierarchy mass problem 
and a realization of electroweak symmetry breaking by 
radiative corrections d]. The typical range of neutralino 
masses in these models extends from a few tens of GeV 
up to - 10 TeV. 

Two notable non-thermal DM candidates are axions 
and superheavy DM (SHDM) particles. Axions were 
proposed as solution to the strong CP problem, but 



the still viable "axion window" includes the possibility 
that axions are the main contribution to the DM abun- 
dance. Depending on the inflationary scenario, axions 
may be non-thermally produced either by the misalign- 
ment mechanism or the decay of axionic topological de- 
fects [i]. 

We will review the status of the SHDM model in the 
next section, before we discuss "superheavy supersym- 
metry" in Sec. HI. Results for the relevant cosmological 
processes of a superheavy neutralino like elastic scatter- 
ing on light fermions, self-scatterings and annihilations 
are presented in Sec. IV. After that, we review various 
possible production mechanisms of SHDM for the special 
case of a superheavy neutralino in Sec. V and summarize 
in Sec.VI. 



II. SUPERHEAVY DARK MATTER 

The first proposal of SHDM in Refs. 0, i] was mo- 
tivated by observations of ultrahigh energy cosmic rays 
(UHECR), which revealed not the expected suppression 
of the energy spectrum due to the interaction of extra- 
galactic protons with cosmic microwave photons. Since 
CDM is gravitationally accumulated in the halo of our 
galaxy, the secondaries produced in decays or annihila- 
tions of SHDM do not suffer energy losses and their en- 
ergy spectrum is characterized by a flat spectrum up to 
the kinematical cutoff mx/2. In this scenario, the mass 
mx of the SHDM particle should exceed 10^^ GeV. Mean- 
while, this original motivation for proposing SHDM has 
disappeared in the light of new UHECR data compatible 
with the expected flux suppression, for more details see 
Ref. 0. 

Superheavy dark matter particles with the density re- 
quired by cosmological observations can be efficiently 
produced at inflation b y m any mechanisms including 
thermal production 0, [8l.ll0l| . The most detailed descrip- 
tion of this process within an inflationary framework is 
given in Ref. plj . 

A variety of different production mechanisms can pro- 
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vide a non-thermal distribution of superheavy particles in 
the expanding universe. Since the energy density of non- 
relativistic particles decreases slower than the one of radi- 
ation, their abundance increases by the factor a{to)/a{t^,) 
with respect to radiation, where a(to) and a(i*) are the 
scale factors of the universe today and at the epoch of 
particle generation, respectively. If particle production 
happens at the earliest relevant time, i.e. during infla- 
tion, this factor can become extremely large, ~ 10^^. 
Not surprisingly, such a small energy fraction can be 
transferred to SHDM particles by many different mech- 
anisms, as thermal production at reheating 0, the 
non-perturbative regime of a broad parametric resonance 
at preheating [13, [13] , and production by topological de- 
fects [3, [11. 

We discuss first the generation of superheavy particles 
by gravitational interactions from vacuum at the end of 
inflation [isl . [l6j . Since this production mechanism relies 
only on the gravitational coupling of the SHDM parti- 
cle, it is unavoidably present in contrast to other, more 
model-dependent generation mechanisms. Neither infla- 
tion is needed for this production, it rather limits the 
gravitational production of the particles. Since this pro- 
duction is caused by the time variation of the Hubble pa- 
rameter H(t), only particles with masses mx ^ H{t) can 
be produced. In inflationary scenarios H{t) < m^, where 
is the mass of the inflaton. It results in the limit 
on the mass of the produced particles, mx < 10^'^ GeV 

The numerical calculations of Ref. p£] for the present 
abundance of fermionic SHDM can be approximated as 

n h'^ ~ / i^x/Hif , mx <^Hi 

^ ~ IQSGeV \ exp(-mx/-ff/) , mx^^Hi ' 

(1) 

where Hj k, 10^^ GeV is the Hubble parameter during 
inflation, and Tr is the reheating temperature. 

Other generation mechanisms can occur additionally 
to gravitational production. If the SHDM particles cou- 
ples directly or through an intermediate particle to the 
inflaton field, the time-dependence of the classical infia- 
ton field induces particle production. SHDM particles 
may be also efficiently produced at preheating [12| . This 
stage, predecessor of reheating, is caused by oscillations 
of the inflaton field relative to the minimum of the po- 
tential after inflation. Such an oscillating fleld can pro- 
duce non-perturbatively in the regime of a broad para- 
metric resonance intermediate bosons which then decay 
to SHDM particles. The mass of the SHDM particles can 
be one or even two orders of magnitude larger than the 
inflaton mass. 

Another mechanism is the so-called instant preheat- 
ing [isj . It works only in specific models, where the mass 
of the intermediate boson x is proportional to the inflaton 
field, — gcj). When the inflaton crosses the potential 
minimum (j) — x particles are massless and they are 
efficiently produced. When increases, m^ increases, 
too, and can reach values up to the Planck mass. 



While these additional production mechanisms can 
increase the abundance of SHDM particles relative to 
Eq. ([T]), entropy production as for instance in thermal 
inflation can reduce their abundance. Therefore there 
exists only a lower, not an upper limit for the mass of 
SHDM, arising from the condition that the SHDM par- 
ticles do not reach chemical equilibrium. 

What are the particle candidates for SHDM? 

The first problem one meets is the particle life-time. 
Superheavy particles are expected to be very short-lived: 
Even gravitational interactions, e.g. described by dimen- 
sion 5 operators suppressed by the Planck mass, result 
in lifetimes much shorter than the age of the universe to- 
Superheavy particles must be thus protected from fast 
decays by a symmetry which is respected even by grav- 
ity. Such symmetries are known: They are discrete gauge 
symmetries. These symmetries can be very weakly bro- 
ken, e.g. by wormhole \f\ or instanton effects to pro- 
vide a sufficiently long lifetime t, t > t^. A systematic 
analysis of broken discrete gauge symmetries is given in 
Ref. For instance, the lifetime of SHDM with mass 
mx ~ 10"'^'^ — 10^^ GeV was found to be in the range 
10^^ — 10^^ yr in the case of the symmetry group Ziq. 

Various models that contain either absolutely stable or 
unstable particles with life-times larger than the age of 
the universe have been discussed [ij, [H, [l^. Most of 
the suggested SHDM candidates belong to a new sector 
that has no tree-level interactions with SM particles. By 
contrast, we study in this work the possibility of having a 
SHDM particle with SM-likc couplings to the weak gauge 
bosons. Since the longitudinal part of gauge bosons cou- 
ples as cx gMx/mz to a particle with mass Mx, weak 
interactions become generically strong for Mx ^ mz 
and thus the perturbative expansion fails. Using partial- 
wave unitarity, Chanowitz, Furman and Hinchliffe (20| 
derived thereby an upper limit of Mx ^ TcV for parti- 
cles coupling with SM strength to the weak gauge bosons. 
An exception to this bound are supersymmetric theories, 
if only mass terms are added that break supersymmetry 
(SUSY) softly [2l|, and in particular the minimal super- 
symmetric extension of the SM (MSSM) [22]. Therefore 
we are led to suggest superheavy supersymmetry, i.e. the 
case where all masses of supersymmetric particles are of 
order 10^^ GeV or larger, as a concrete model for SHDM 
with SM weak interactions. For definiteness, we choose 
the LSP as the lightest neutralino but we note that other 
possibilities as a sneutrino are also viable. 



III. SUPERHEAVY SUPERSYMMETRY 

Our main motivation for the introduction of super- 
heavy supersymmetry is the search for a particle can- 
didate for superheavy dark matter. As discussed in 
the previous section, the particle candidates found so 
far [13, [H, [3 particle sectors, such as e.g. 

the hidden sectors of supergravity or string models. In 
this section we discuss a more natural and more familiar 
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candidate, the lightest supersymmetric particle, in the 
case that all supersymmetric particles are superheavy. 
The longevity of this dark matter candidate is provided 
by a Z2 discrete gauge symmetry (R-parity) and its pro- 
duction is guaranteed by gravitational interactions in any 
standard inflationary scenario. The scale of SUSY break- 
ing may be determined by LHC experiments, and our 
model can be soon falsified by the discovery of low-scale 
SUSY at LHC. 

Independent of the scale of symmetry breaking, super- 
symmetry remains an inevitable feature of any theory 
that wants to unify internal gauge symmetries such as 
SU(5) or SO(IO) with the symmetry group of Minkowski 
space-time, the Poincare group. Promoted to a local 
symmetry, gravity and gauge interactions are on a sim- 
ilar footing, with the gravitino as gauge field of gravity. 
Another motivation to consider (superheavy) SUSY is 
that it is a generic ingredient of consistent string theo- 
ries. Thus various esthetical reasons suggest that SUSY 
may be realized in Nature. 

Below we discuss the status of superheavy supersym- 
metry in comparison with low-scale symmetry breaking. 

There are three pieces of evidence pointing towards 
"low-scale SUSY" , i.e. a mass scale Msusy of the super- 
symmetric partners of the SM particles below or around 
1 TeV. First, the unification of coupling constants fails 
in the SM, while the three couplings meet in the MSSM 
assuming AfsusY ^ 1 TeV [23J. Second, the fine-tuning 
problem of the SM Higgs is remedied in the MSSM, only 
if the mass splitting between the SM particles and their 
SUSY partners is small enough. Third, the SM contains 
no suitable DM candidate, while the lightest supersym- 
metric particle (LSP) of the MSSM, assuming that R 
parity is conserved, is a promising CDM candidate. As- 
suming further that the LSP is a thermal relic requires 
again that at least part of the SUSY mass spectrum is 
below or close to the TeV scale. 

These attractive properties of low-scale SUSY are over- 
shadowed by several less appealing features: Low-scale 
SUSY predicts generically excessive fiavor and CP vio- 
lation as well as proton decay through dimension-5 op- 
erators that is close to or exceeds observational bounds. 
Moreover, the non-observation of SUSY particles with 
masses around 100 CeV re-introduces a "small" fine- 



tuning problem [24| : For instance, electroweak symmetry 
breaking requires that 



m? tan^ 



/3 



tan" /3 — 1 



(2) 



where mi are the usual mass parameters of the Higgs po- 
tential that depend implicitly on the soft SUSY breaking 
masses. For <C mf^fj?, a certain amount of cancel- 
lation between the terms on the RHS is required. Simi- 
larly, an upper limit on neutralino mass mass of 200 GeV 
arises, if one limits accidental fine tuning to the level of 
1% [2^. Moreover, most of the so-called "bulk region" 
in which the lightest neutralino has naturally the correct 
DM abundance is meanwhile excluded I2al. 



Recently, split SUSY has been proposed as a model 
avoiding the problems of low-scale SUSY while keeping 
gauge coupling unification [13] ■ In this model, the mass 
spectrum of SUSY particles is separated in two parts: 
Gauginos and gluinos are kept at the TeV scale provid- 
ing with the lightest neutralino a suitable thermal DM 
candidate, while all scalars additional to the SM Higgs 
are heavy, with masses possibly close to the CUT scale. 
Motivated by the cosmological constant problem and the 
landscape picture [11] suggested by string theory, the nat- 
uralness principle is given up, keeping as guiding prin- 
ciples only experimental observations: The existence of 
DM, and the hint for CUT from gauge coupling unifica- 
tion. 

In this work, we go one step further by abandoning 
also for the gaugino and gluino masses the weak scale. 
This becomes possible because we assume that the LSP 
is produced non-thermally at the end of inflation. As a 
result, the mass of the LSP should be generically above 
- 3 X 10" CeV. The mass of the gluinos and of the SUSY 
scalars could be either close or, in a similar but not as 
extreme set-up as in split SUSY, much larger. In such a 
set-up, alternative approaches as modular [29| or confor- 
mal [30| invariance may provide a solution to the hierar- 
chy problem and to gauge coupling unification. 

Experimental data from LHC will soon decide if the 
supersymmetric particles are at least partly close to the 
weak scale. If this is not the case, then both the MSSM 
and split supersymmetry are disfavoured. Our proposal 
that superheavy LSPs (SHLSP) are the DM particles 
may be then an interesting alternative connecting SUSY 
to the physical world. The prospects to detect DM in the 
form of stable superheavy neutralinos will be discussed 
in a subsequent work [sij . 



A. Neutralino as LSP 

We assume throughout that the lightest neutralino x = 
Xi is the lightest of the supersymmetric particles. The 
neutralino mass-matrix in the {B, H^, H2) basis 
is given by [s^] 



Ml 


-cpmzsw 
mzSwSp 





M2 
cwcpmz 
-cwmzSjS 



-cpmzsw 
cwcpmz 




mzswsp 
-cw'm.zSfi 





(3) 

with ~ sin/3, C/3 — cos/3 where tan/3 — Vi/v2 is the 
ratio of the two Higgs vev, sw = siwdw, cw — cos-dw 
with i!)\Y as Weinberg angle, and fi as the Higgs mixing 
parameter. 

The consequences of the limit Mi,M2, ^ mz for 
the neutralino has been already extensively discussed for 
a neutralino as thermal relic [33] . Neglecting the terms of 
order mz , the four neutralino mass eigenstates become a 
pure bino, wino, and the symmetric and anti-symmetric 
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combination of the two higgsinos, 

{B, M^", (i/o + H',)/V2, (i/o - i/°)/V2} (4) 
with masses 

{^h,M2,-^l,^^}. (5) 

In order to decide which of the two higgsino combinations 
is the lightest, one has to include corrections of second 
order in mz- Then one obtains as mass eigenvalues 



Ml 



Mo 



si, (Aisin(2/3) + Afi)m| 



Mf 



(sin(2/3) - 1) (a* + c^Mi + s^Mj) m 



2(m + Mi)(m + M2) 
(sin(2/3) + - c^Mi - g^M2)m^ 
2(Ai-Mi)(^-M2) 



}, (6) 



if the masses are not degenerate. Depending on the 
sign of /X, the symmetric {fj, < 0) or the anti-symmetric 
(/U > 0) combination of the higgsino is the LSP for 
^ Mi,M2. To be definite, we shall choose always 
/i > in the following. 

If the mass difference between the two lighter of the 
three mass parameters ^, Mi and M2 are small com- 
pared to TO^, a large mixing between the gaugino and 
the higgsino remains even in the limit |/i|,Af 3> mz- 
This case of a (partially) degeneration between the neu- 
tralino mass parameters /i. Mi and M2 was dubbed well- 
tempered neutralino and discussed in Ref. [s^l- 



B. Unitarity for superheavy particles 

We illustrate with two explicit examples how super- 
heavy SUSY avoids that longitudinal gauge bosons be- 
come strongly coupled to neutralinos for ^ mz- The 
first case is the annihilation of neutralinos into fermion 
pairs in the limit of zero relative velocity v. Then 
the amplitude consists of sfermion, Z and A exchange, 
A4 = Aiz + -Ma + -Msf- An inspection of the annihi- 
lation cross section given e.g. in Ref. [4| shows that the 
longitudinal component contributes the term 



1l 

IGtt 



g 



Oi^?T. 



2 



(7) 



to the annihilation cross section into fermions with mass 
m/, isospin Ta, couphng O^f , and (3j = l — m^^/m'^. 
Assuming that there are no cancellations between A4z 
and Ma+MsJ as well as that the factor Oii is of 0(1), 
the annihilation cross section would be independent from 
and thus violate perturbative unitarity for 3> 

mz- 



The resolution of this apparent problem lies partly in 
the specific form of the neutralino coupling and mass 
matrix. The interaction between Z^ and the neutralinos 
are given in the unitary gauge by the Lagrangian Q 



-Zx'>X° 



9 



2cw 



ZMno'rlPL + 0'2Pr)x'L], (8) 



where n — m — 1 yields the interaction between and 



the LSP. Moreover, expressing O-^ ' by the neutralino 
mixing matrix elements gives 



0,1 = ~Oj: = U-NsnN;,^ + N.^NU . (9) 

Neglecting as always below possible CP violation, the 
coupling becomes simply cx (iVfa — for m = n = 1. 
For a completely (anti-) symmetric higgsino, this cou- 
pling vanishes because of jA^is] = |iVi4|, while a bino 
and wino LSP have jiVial = |A^i4| = in the limit 
mz 0. The leading contribution is thus suppressed 
by (m^/Msusv)^ and given by 



s(2/3)s 



2(m"-M|) 
cos(2/3)m|(c' 



if Ml < A/2,M, 

if Af2 < Mi,n, 



,Mi+s^M2-fi) 



4M(M-Afi)(M-Af2) 



if Ai< A//i,A//2. 

(10) 

Thus this coupling vanishes in the limit m^/mz 00, 
because the longitudinal components of the gauge bosons 
couple only to the deviation from a completely (anti-) 
symmetric mixing of the higgsino components. Although 
the approximations (jlOp are valid only for non-degenerate 
masses, this conclusion holds also for degenerate masses, 
because lA^ia] « |A^i4| remains valid. Taking into ac- 
count this suppression factor from the neutralino mixing 
matrix, already the single term from Z^ exchange in the 
annihilation cross section is consistent with the unitary 

bound, (Tann OC 

Another way to understand how the apparently dan- 
gerous terms m^jmz disappear in physical quantities is 
to compare the coupling of the Z and its goldstone boson 
Gz in different gauges. For simplicity, we use for this 
comparison the unitary and the i?j-gauge restricted to 
^ = 1 . Then the Z propagator becomes purely transver- 
sal in the i?^-gauge, and the interactions of neutralinos 
with the goldstone Gz have to agree with those with the 
longitudinal part of the Z boson in unitary gauge. 

We consider as example the annihilation of neutralinos 
into a Z and the lighter CP even higgs boson in the limit 
of vanishing relative velocity v. Then the longitudinal 
part Zj gives for x(p)-|-x(p') Z{k) + h{k') annihilation 
at rest |3a| 



M{xX^Zh)z, 



cos^ vw rnz q 



(11) 



where q = p + p' denotes the momentum of the virtual 
Z]^ and e the polarization vector of the real Z boson. 
Although the neutral Goldstone boson is part of the Z 
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boson in the unitary gauge, its coupling to neutralinos 
differs and is given by 



G 

XnXr, 



2cosi? 



W 



(12) 



with 



O 



G 

nm 



(13) 

For the lightest neutralino annihilation, n — m — 1, and 
since the coupling is imaginary, the Goldstone only cou- 
ples to the axial part of the neutralino. The amplitude 
for the Goldstone exchange diagram is thus 



M{xx^Zh)c 



q-ejk) 



(14) 



Comparing (jlip and (jl4p . gauge independence requires 
as relation between the couplings 



'-^11 



mz 



(15) 



or expressed in terms of the neutralino mixings and 
masses, 



IV. RELEVANT PROCESSES AND CROSS 
SECTIONS 



Elastic scattering on fermions and the energy 
relaxation time 



Kinetic equilibrium of neutralinos in the late universe 
may be reached by scattering on light fermions like neu- 
trinos and electrons. In the rest frame of the neutralino, 
the Mandelstam variables become 



t = -2tj2(l - cosd) , 



(19) 



where is the mass of the lightest neutralino, lo is the 
initial energy of the lepton and d is the scattering angle. 
We consider here only the case of a broken electroweak 
symmetry, i.e. the case of temperatures T below the weak 
scale, when the following hierarchy holds 



UJ <^ mz ^ , 



(20) 



The assumption mz <C m^^ leads also to several simpli- 
fications in the Higgs sector of the MSSM that we shall 
employ below. Additionally, we require that the neu- 
tralino mass parameters are not too degenerate. 



mz 



'{cwNi2 - swNn){spNi4 + cpNis) . (16) 



The authors of Ref. [35| noted that this relation can be 
derived directly from the definition of the neutralino mix- 
ing matrix N, 



(17) 



since the identity 



(Wri 



m,„) (iVFiV-^)„,„ = {N (MP + PM) 7V-i)_ 

holds for any matrix P. Choosing for P the isospin op- 
erator that flips the first higgsino sign compared to the 
second one, P = diag(0, — 0-3), reproduces a generaliza- 
tion of Eq. (fTB]) . It is the special structure of the neu- 
tralino mass matrix that makes PM -\- MP off-diagonal, 
and thereby leads to the vanishing of the left-hand side 
of mi) with mz- 

Finally, superheavy particles in a theory with chiral 
particles may lead to radiative effects that do not van- 
ish for 00. The authors of Ref. [l^l discussed in 
a series of works, if the SM can be viewed as the low- 
energy limit of the MSSM in the sense of the Appelquist- 
Carazzone theorem [3^. They showed that all virtual 
effects of the SUSY particles are either suppressed by 
inverse powers of their mass or can_be absorbed in the 
renormalization of SM parameters 



4l|. In conclusion. 



superheavy SUSY particles do neither lead to a violation 
of perturbative unitarity or to non-decoupling effects in 
virtual corrections. 



|/i-Mi|, \n-M2l IM2-M1I >mz 



(21) 



Wc consider explicitly the case where the lightest neu- 
tralino is a bino or a higgsino and scatters on a neutrino. 
The case of a wino is almost identical to the one of the 
bino. 

The Feynman amplitudes of the process x4-J^e ~> X+^^e 
consists of three contributions: Sneutrino exchange in the 
s and u channel, and t channel exchange of higgses and 
the Z, \M\'^ = \Ms — Mu + A^tp. Since the neutralino 
is a Majorana particle, the amplitudes Ms and Mu can 
be obtained by interchanging the neutralino in the initial 
and final states and thus they should be subtracted. 



1. The bino as the LSP 

Using the approximations explained above, we obtain 
as the leading contribution to the total spin-av erag ed 
squared Feynman amplitude in the case of a bino 



\Mu\^ = \Ms 



e^Mfu;^ 



(22) 



2 e4MV(3-cos(.9))cos^(2/?) 

1-^*1 ~r~^ — r::TT2 > i^^^i 



2c4, (^2 „ 



2MMsM:)^-q^l^f^, (24) 
ctv (M? - Ml) 
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2ReiMsM;) = -2ReiMtMl) = 

cos(2/3)Mi2w2(3 - cos(z?)) 
2c4^(M2-M?)(Af2_^2) 



(25) 



Here, we used neutrinos as scattering target and denoted 
by Mp the sneutrino mass. 

The e nerg y relaxation time can be calculated as (see 
e.g. Ref. [11]) 



1 

Trcl 



iVcff 



(26) 



where Ek ~ {3/2)T is the mean kinetic energy of the 
neutralinos, dp the ncutralino momentum obtained in one 
scattering, 



{Spf = 2uj^il-cos{d)) 



(27) 



and the number density of relativistic fcrmions with one 
polarization and energy w is 



1 



no 



27r2 



1 



"/T + 1 27r2 



cj^ e 



(28) 



Finally, the factor A'cff counts the number of relevant rel- 
ativistic degrees of freedom, weighted with the relative 
size of their cross-section compared to a neutrino. Com- 
bining the different contributions and performing the in- 
tegrals gives 



Trcl 



25iVoffe4r6 



cos(2/3)Af| -f /i2 - 2c^M2 



(29) 



2. The hiqqsino as the LSP 



In the case of a higgsino as the LSP the contributions 
to the total squared Feynman amplitude \M\^ are given 
by 



\Mu\' = 



(30) 



(31) 



\Mt 



2e^up-c\i^ (3 ~ cosi^) (Mic^ + A^s^r - n)' 



4^ {fx-A'hfip-M.f 



2Re{MsM*t) - 2e^ (c2/3 + Sip/2) (3 - cos^) 
4^(M-Mi)^(M-M2)3(Ai2_M?) 



(32) 



(33) 



2Re{MsM*J = -46* (c^ + spf sm''{d/2) 

s\^{^i-^h)\^Ji-M^)\Ml-^i^f 



(34) 



2Re{MtMl) - -2e^ {c^p + Sip/2) (3 - cosi?) 



(35) 



Analogously to the case of the bino, the energy relaxation 
time for a Higgsino as lightest neutralinos follows as 



Trcl = 



ttV^ {^l-^hy{^l-M2ys■ 



2W 



(36) 



B. Elastic neutralino-neutralino scattering 

We use in this subsection again the assumptions (|20p . 
but denote now with lu the kinetic energy of the collid- 
ing neutralinos in their center of mass frame. Then the 
Mandelstam variables are 

5 = 4(1/2-1-^2), t = -2^2(1 -cost?). (37) 



1. The Bino as the LSP 

Neutralino-neutralino scattering can occur through s 
and t channel exchange of the Z and the three neutral 
Higgs bosons. We shall see that it is sufficient to consider 
only the squared matrix elements. They are given in the 
unitary gauge by 



\Mz-a,nn\ = 



2 9e^ cos^{2f3)m%M^ tan^ jdw) 
~ 2{„^-Mff 



2{^,^~Mlf 
e^cj^ml (^ sin(2/3) -f Mif tan^ [dw) 



(38) 
(39) 



2Af4 (At2 - Aff; 



(40) 



\M 



h-cxch 



Be^mlM^ (^ sin(2/3) + Mi)* tan* {dw) 



Mt (^2 - A/2; 



\Mh- 



cxch I 



,2 _ 8eVcos*(2/3)m|A/ftan*(t9vF) 



\Mh-. 



A/4 (;,2 _ Ml)' 

jeV^cj* cos*(2/3)m| tan* {dw) 

{Ml - AMlf {p^ - A/2)* 
Be*m|Mi* {p + sin(2/3)A/i)* tan" (t9n/) 



(41) 
(42) 



(43) 



(44) 
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and 



\M 



A"-GXChl 



(45) 



u;^(cos(27?) + 7)m| (^ + sin(2/3)Mi)^ 
M4„ (^2 _ Mlf ■ 

We note first that we can neglect the squared amphtudes 
proportional to w'*. Because of the hierarchy in the Higgs 
masses, 



0{Mh) = 0{mz) < 0{Mh) = 0{Mao) 



(46) 



the /i-exchange channel (|¥T|) that is of order 0{Mg^gY) 
compared to other channels of 0{m'^ /Mg^gy) domi- 
nates the self-scattering of superheavy neutralinos. With 



'X"X"^X"X"l — l-'^/i-cxch 

neutralino-neutralino scattering follows as 
_ e^mlMf (^sin(2/3) + Mi)^tan4 



the total cross section of 



(47) 



2. The Higgsino as the LSP 

Analogously to the bino case, the leading contribution 
to higgsino-higgsino scattering is given by the exchange 
of the light, SM-like Higgs h. With 



\Mh-. 



lJL^m% (^1 cos^ {'dw) - + sin^ {dw) M2) 



(48) 



the total cross section of neutralino-neutralino scattering 
follows as 



256^ct^s4,M4 (^ - Ml)" (a* - M2)" 



(49) 



C. Annihilations 

The annihilations of neutralinos have been studied 
in great detail. Annihilations of superheavy neutrali- 
nos are, in the bino case, dominated by the channels 
ZH, hA, AH, W^H^, since all other channels are sup- 
pressed by powers of m^/MsuSY- In the higgsino case 
all bosonic channels contribute at leading order except 
annihilation into + A^ and h + H. Annihilation into 
fermions are always suppressed. 

The fermionic annihilation channels do not give leading 
order contributions in any case. 



The bino as the LSP 



these dominant channels are given by 

{-A^lMf + A + B + fi^Mls2f^y 



(50) 



8cty (m^ - MfY [Mi - 4Mfr (-2^2 - 2M/ + Mf,) 



with 



A = {Ml-2Ml)s2i3Ml 
B = 2fi (2^2 ^ ^,j2 _ ^j2^-) 



(51) 
(52) 



\M^o^o^aoh\' = (Ml - M],y /{ScUh) x (53) 
{-&S2f3Mf - SfiMf + C + D + Ef 



ip^ - Mff {Ml - AMlf (-2^2 _ 2Af2 + Ml + Mjjf ' 



with 



C 
D 
E 



2 {Ml + Ml) S2pMl 



= 4^ (-2^^ 



Mi 



Mfj) Ml 



cpMl + Ml + Ml) sp 



I , . ,2 (sin(2/3)Ml-f4MMi) 

'''' 8c4^ (Ml- 2^2 _2Mf)' 



and 



(54) 
(55) 
(56) 

(57) 
(58) 



IM^o^o^w^htI'' (M^± -4M2)'/8c4^ > 
(-4^M3 + F + G + fi'^Mls2p)'^ 
(M^± - 2^2 - 2M2)' (^2 _ {Ml - 4M2; 



where 

F = {Ml-2Mjj±)s2pMl (59) 
6 = 2/1 {-Mjj± + 2/^2 + Ml) Ml . (60) 

The annihilation cross section of the various channels fol- 
lows then as 



327rm2 



(61) 



2. The higgsino as the LSP 

Assume that {M^-^, M^^, M~^^, M^^} correspond to 
— /i. Ml, M2}. Then the squared matrix-elements of 
these dominant channels are given by 



\M 



2e^ 



sin {2i} 



w 



\4 ' 



(62) 



Assume that {M^-^, M^^, M^^, M^^} correspond to 
{Ml, M2, /i, — /i}. Then the squared matrix elements of 



4sin(7?iv)* ' 



(63) 
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\M^o-^a^Z«h\^ 



•2/3 



if,^ + Mir {^^^ + Mir 



X (64) 



(c/3 + S/3) 



2W 



(ifiMislrA - 2clvMlB + 4/^c^M2C + Z?)^ ,(65) 



where 



A = (-2^2 _ 2M| + M^) (66) 

B = (Ml + 4/iAf2) (67) 

C = (Ml -2/^2) (68) 

1} = M\~2^?M\~2MIswMa. (69) 

The annihilation into , ^ are dominated by the 
chargino exchange and the A^ annihilation channels. In 
order to obtain relatively compact expressions, we set 
Mj^o = Mh± and assume that ficos{f3) + sin(/3)M2 ^ 0. 
As a result we obtain 



{s2pMfj^ + 4^M2) 
8 (-2^2 _ 2M2 + M|±) 



.(70) 



To simplify the squared matrix elements both for anni- 
hilations into + iJ and A" + H, we set M^o ^ Mh- 
The results are 



icf3 - 8/3) ' 



'r'r'-^^-n , ^ ^_2^2 _ + M^) 2 
(M^ - 2^2 ^f^ + S - F - G) ^ (71) 



with 



E = 4AiMi (2^2 _^ 2Af2 - M|r) (72) 
F = 2c2^Mf (M|f - 4AiM2) (73) 
G = 2M2{2fi{Mfj-2^j,^)clr + M2Mljsiv) (74) 



H = 2s. 



2W 



^2^2 _ 2M2 + Mfj) 



(75) 



and 



2eV'c? 



|Mx»X«^AO«|^ - (^2+M2-?f^)2/ ^ 

(Mf M2C^ + Ah (Ai' - M^) + J) ' (76) 



V. COSMOLOGICAL PRODUCTION 

All the mechanisms for the production of SHDM par- 
ticles described in section |lT] apply also to SHLSPs. In 
particular, the abundance of neutralinos due to gravi- 
tational production is given by Eq. ([1]). In most works 
on SHDM particles the dependence on and the two- 
fold degeneracy in Eq. ([T]) was fixed by choosing first 
for Tji the highest value allowed by the gravitino prob- 
lem, Tfj — 10^ GeV. Then the larger of the two possible 
masses was selected, ~ 3 x 10^'^ GeV, so that secon- 
daries of SHDM decays could explain the cosmic rays of 
the highest energies, E > 10^^ eV. 

Both constraints can be relaxed in our scenario: Super- 
heavy gravitinos decay before big-bang nucleosynthesis 
and we do not insist that SHDM decays explain UHE- 
CRs with energies, E > 102°eV. As a result, the only 
constraint additionally to Eq. ([T]) is the requirement that 
neutralinos do not thcrmalize, <§; m^/30. Thus the 
choice Tb. = 10^° GeV and = 3 x 10" GeV gives 
the smallest possible neutralino mass. Finally we note 
that decays of heavier SUSY particles with mass rh, that 
are more effectively produced in the regime m <^ Hj 
than neutralinos, can change the relation ([T|) but not the 
smallest possible neutralino mass. 



VI. SUMMARY 

We have suggested the lightest supersymmetric parti- 
cle as a well-suited candidate for superheavy dark matter. 
The requirement that SHDM is produced with the cor- 
rect abundance by gravitational interactions at the end 
of inflation leads to a lower mass limit of 3 x 10^^ GeV 
for the masses of all SUSY particles. Since weak interac- 
tions remain perturbative despite the large mass hierar- 
chy, 3> mz, and the mass scales are approximately 
fixed, we could reliably calculate the relevant processes 
for the special case of a superheavy neutralino. 

Our proposal can be falsified in the near future by the 
discovery of low-scale/split SUSY at the LHC. If this 
is not the case, then SHLSPs as DM particles may be 
the unique opportunity to connect SUSY to the physi- 
cal world. The observable consequences of meta-stable 
SHDM were discussed already in detail in Refs. ^4^, 
while the prospects to detect DM in the form of sta- 
ble superheavy neutralinos despite of their small number 
density and annihilation cross section will be discussed 
in a subsequent work [3l| . 



with 



'2W 



/ = {^? + ^q-Ml) 



(77) 
(78) 



The annihilation channels Z° -I- A'^ and h + H do not give 
an 0(1) contribution. 
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